Entanglement probe of two-impurity Kondo physics in a spin chain 
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We propose that real-space properties of the two-impurity Kondo model can be obtained from an effective 
spin model where two single-impurity Kondo spin chains are joined via an RKKY interaction between the 
two impurity spins. We then use a DMRG approach, valid in all ranges of parameters, to study its features 
using two complementary quantum-entanglement measures, the negativity and the von Neumann entropy. This 
non-perturbative approach enables us to uncover the precise dependence of the spatial extent £x of the Kondo 
screening cloud with the Kondo and RKKY couplings. Our results reveal an exponential suppression of the 
Kondo temperature Tk ~ l/£ k with the size of the effective impurity spin in the limit of large ferromagnetic 
RKKY coupling, a striking display of "Kondo resonance narrowing" in the two-impurity Kondo model. We also 
show how the antiferromagnetic RKKY interaction produces an effective decoupling of the impurities from the 
bulk already for intermediate strengths of this interaction, and, furthermore, exhibit how the non-Fermi liquid 
quantum critical point is signaled in the quantum entanglement between various parts of the system. 

PACS numbers: 71.10.Hf, 75.10.Pq, 75.20.Hr, 75.30.Hx 



Introduction.- The theory of quantum impurities underpins 
much of the current understanding of correlated electrons. A 
case in point is the two-impurity Kondo model (TIKM) U], 
with bearing on heavy fermion physics J2l, correlation effects 
in nanostructures J^, spin-based quantum computing fl,|5t], 
and more. The model describes two localized spin- 1/2 im- 
purities in an electron gas, coupled by the Ruderman-Kittel- 
Kasuya-Yosida (RKKY) interaction ~ Jj via their spin ex- 
change with the electrons. In addition to the RKKY cou- 
pling the model exhibits a second energy scale, the Kondo 
temperature Tk, below which the electrons may screen the 
impurity spins. For strong ferromagnetic RKKY interaction, 
| Jj \ Tk, the impurities form a spin-1 state which does get 
screened, in exact analogy with the spin-1 two-channel Kondo 
effect. In contrast, for strong antiferromagnetic RKKY inter- 
action, J i Tk, the impurity spins form a singlet state, 
killing off the Kondo effect. In the presence of a special 
electron-hole symmetry ffj, or with the impurities coupled 
to separate electron reservoirs J3l, the nonuniversal crossover 
between the two regimes sharpens into a quantum phase tran- 
sition (QPT) with a non-Fermi liquid quantum critical region. 

While much is known about the model, various real-space 
properties are yet to be uncovered. A central concept is the 
Kondo screening cloud, inferred from the appearance of a 
characteristic length £x = Hvp/ksTK associated with the 
energy scale Tk, with vp the Fermi velocity HI]. The Kondo 
cloud is invoked to explain how physical quantities at a dis- 
tance r turn into scaling functions of t/£k, rather than de- 
pending on r and £k separately 181]. However, its nature, 
structure, and experimental confirmation has remained con- 
troversial, even for the single-impurity Kondo model, moti- 
vating several recent attempts to determine the Kondo length 
£k accurately JHQll. The case of the TIKM is further com- 



pounded since the Kondo regime changes its character as one 
tunes through zero RKKY coupling (where the model splits 
into two single-impurity Kondo models) to strong ferromag- 
netic RKKY coupling (where the physics is instead captured 
by the spin-1 two-channel Kondo model). What is the signa- 
ture of this crossover? How does it show up in real space? 

In this Letter we address these questions by exploiting and 
making precise the picture that a screening cloud is built from 
those electron states that are entangled with the impurities 
Jit]. This allows us to nonperturbatively uncover: (i) the quan- 
tum phase transition between Kondo and RKKY regimes; (ii) 
the effective impurity-bulk decoupling for antiferromagnetic 
RKKY coupling; (iii) the true spatial extent of the Kondo 
cloud; and (iv) the effect of "Kondo resonance narrowing" for 
large ferromagnetic RKKY interaction. 

The spin model.- Exploiting the effective one- 
dimensionality of the TIKM iflUl . we introduce its spin 
chain emulation by couplin g th e impurities of two single- 
impurity Kondo spin chains ll 211 by an RKKY interaction of 
strength Jj, see Fig. [Ha). The speed-up of numerics achieved 
by working with a "spin-only" version of the TIKM is 
significant, and enables us to extract entanglement properties 
via a high-precision Density Matrix Renormalization Group 
(DMRG) approach. We thus consider the spin Hamiltonian 
H = Y,k=LR H k + H T where 
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Here k = L,R labels the left and right chains, with cr\ the 
vector of Pauli matrices at site i in chain k, and with Ji ( J-f) 
nearest- (next-nearest-) neighbor couplings. Taking J\ = 1, 



2 



(a) LeftBulk 



Right Bulk 



N N-1 



2 3 



N-1 N 



\/\/\ / \| }--V-f- j \/\/\7\ 



(b) 



\/\/ 


\/\ 


I V-f 


\/\ 


/\/\ 


(c) 


M-2 








\/\/ 




T~ ^ 


\/\ 


/\/\ 



«0 ^ 



Traced out spins 



Traced out spins 



FIG. 1 : (Color online) (a) The composite system of two Kondo spin 
chains, each formed by an impurity and a bulk, and coupled together 
by an RKKY interaction ~ Ji to represent a "spin-only" TIKM. 
(b) The symmetric central block which contains M spins, including 
impurities, (c) The effective model in which the two impurities are 
decoupled and instead the bulks are interacting through Hf , given 
in Eq. {4j. (d) Two identical blocks of spins are traced out from the 
bulks and entanglement between the two impurities and the rest of 
the system is computed to find the Kondo length. 



J 2 must be tuned to a critical value Jf = 0.2412 in order 
for Hl and Hr to each faithfully represent the spin sector 
of a single-impurity Kondo model: when J 2 > Jf the spin 
chain enters a gapped dimerized phase, whereas for J 2 < Jf 
a marginal coupling (in the sense of RG) produces logarithmic 
corrections which pollute numerical data Il2lll3ll . The posi- 
tive parameters Jf and J' R play the role of antiferromagnetic 
Kondo couplings. Hi, finally, is the RKKY interaction be- 
tween the impurity spins erf and erf, with Jj allowed to take 
any positive or negative value. In the DMRG calculations we 
take Nl = Nr = N/2 and assume that J' L = J' R = J'- 



Structure of entanglement- For Jj = 0, the system de- 
couples into two independent ("spin-only") single-impurity 
Kondo models where each impurity is screened by its own 
cloud. For nonzero Jr, the ground state attains a more com- 
plex structure. To perform a diagnostic we shall use two 
well known entanglement measures: negativity 11411 and von 
Neumann entropy. Recall that for a bipartite system AB, 
negativity is defined as N(pab) = Yli \ a i\ ~ 1> wnere a i 
denote the eigenvalues of the partial transpose of the den- 
sity matrix pab with respect to one of the two subsets, i.e. 
A or B. The von Neumann entropy, in turn, is given by 
S(p A ) = -TrpA log 2 PA, where p A = TibPab is the re- 
duced density matrix of the subset A when system B is traced 
out from pab- While S(pa) quantifies entanglement only 
when the total density matrix pab is pure, there is no such 
restriction for the negativity. 

Let us first recall that RG studies of the electron-hole sym- 
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FIG. 2: (Color online) (a) Entanglement between two impurities vs. 
Ji for N L = N R = 200. (b) The coupling Jf vs. 1/J'. 



metric TIKM have identified the fixed-point Hamiltonian for 
small positive Jj as that for two independent Kondo impu- 
rities ||2[] . However, unlike for independent spins the ground 
state expectation value (erf ■ erf) is nonzero 11511 . For Ji > 
2.2Tk, no Kondo effect occurs. Still, Kondo correlations per- 
sist in this regime, with the two impurity spins locking into 
a singlet only for very large values of Jj. This picture holds 
also without electron-hole symmetry when the impurity spins 
couple to separate electron reservoirs with no charge trans- 
fer |7|]. As this is the case modeled by our Hamiltonian in 
Eq. (Q]), we expect a quantum phase transition (QPT) at some 



value J i = Jf with Jf scaling with Tk 



,-a/J' 



a being 



a positive constant M15H . To detect the QPT, we compute the 
negativity between the impurities as a function of Jj. One 
may take advantage of the SU(2) symmetry of the system and 
write the reduced density matrix p 1 1r of the two impurities 
as a Werner state, 



Ps\S)(S\ 



l-p. 



(2) 



i=0,±l 



where \S) is the singlet state, |Tj) (i = 0, ±1) are triplets and 
P s is the singlet fraction which varies with Jj and J'. The 
negativity for a Werner state coincides with its concurrence 
IU6I1 and can be obtained as A/"(1l, 1r) = max{0, 2P S — 1}. 
The numerical results are depicted in Fig. |2ja) where the 
entanglement rises from zero (Kondo regime) at a point Jf 
(which depends on J') and eventually saturates to unity (local 
RKKY singlet). Moreover, as seen in Fig. |2jb), the transition 
point Jf indeed scales exponentially with 1 / J', in agreement 
with the RG picture from Ref . |fl5Tl . The small finite size cor- 
rection for small J', captured in Fig.[2|b), reflects the fact that 
in this parameter regime the extent of Kondo screening cloud 
approaches the system size. 

By decreasing Jj, the singlet fraction P s decreases mono- 
tonically and approaches zero in the limit of large negative Jj. 
It follows that in this limit the two impurity spins effectively 
behave as a single spin- 1 entity, with the two bulks serving as 
two screening channels. To see this effect one may calculate 
the von Neumann entropy of p l , 



-P 5 log 2 (P s )-(l-P 5 )log 2 (i-^). (3) 
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As the ground state of the whole system is a pure state, this 
quantifies the entanglement between the two impurities and 
the rest of the system. In Fig. [3ja) we plot S(p 1 1r ) as a 
function of Jj which clearly implies the limiting behavior 
~ log 2 (3) (corresponding to P s = 0) for Jj <C 0. It is in- 
teresting to contrast the Kondo screening behavior for Jj < 
with that for Jj > 0. For Jj = each impurity is maxi- 
mally entangled with its neighboring chain independent of the 
value of J'. By additivity of the von Neumann entropy it fol- 
lows that S(p 1 ! ) = 2, as seen also in the DMRG data in 
Fig. Oa). Returning to Eq. note that P s determines the 
effective impurity spin as a function of Jf. For P s = (i.e. 
J i <C 0) the two impurities behave like a single spin-1 object, 
while for P s = 1/4 (i.e. Jj = 0) the two impurities are de- 
coupled and each, carrying spin 1/2, gets screened by its own 
cloud. In the limit Jj > Owe have that P s = 1, and the two 
impurities form a singlet. 

It is also instructive to study the entanglement between dif- 
ferent constituents of the system. In Fig. [5Jb) we display the 
negativity 7V(1l, -Bl) between the left impurity and the left 
bulk (by symmetry we have N(1 l ,Bl) = N'O-r, Br)). As 
\Ji\ increases, N{1l,Bl) drops rapidly for J/ > as one 
tunes through the QPT where the Kondo screening becomes 
feeble. In contrast, for J/ < the decrease of A/"(1l, -Bl) 
is slower and approaches a finite value in the limit where the 
impurity states form a spin-1 state. However, since the left 
(as well as the right) impurity is now less screened by its own 
bulk, entanglement monogamy Il7tl implies that it is entan- 
gled also with the opposite bulk, as revealed by Fig.fJIc). To 
display the entanglement between the two bulks for finite J] 
we plot the negativity JV(Bl, Br) between the left and right 
chains in Fig.[3jd), having traced out the impurity states. As 
Fig. [3jd) shows, N(Bl,Br) is no longer bounded by unity, 
reflecting the fact that the bulks contains many spins. Further- 
more, due to entanglement monogamy, JV(Bl, Br) is larger 
for Jj > for which the two impurities tend to decouple 
from the rest by forming a singlet, in comparison to Jj < 
for which the two impurities are entangled with the bulks, thus 
reducing their ability to get entangled with each other. 

Quantum phase transition.- To corroborate that Jf is a 
quantum critical point, we plot the first derivative of the nega- 
tivities Bl), AT(li, Br) and M{B L , B R ) inFig.gta)- 
(c) for two cases, i.e. N = 240 and N = 400. The cusps 
at Jf, which become increasingly sharper for larger N, are 
finite-size precursors of a divergence in the thermodynamic 
limit, a hallmark of a second-order QPT IU8I1 . The data in 
Figs. |4|a)-(c), together with Fig. |2|b), provide a highly non- 
trivial check that our spin chain model is a faithful ("spin- 
only") emulation of the TIKM. 

Effective decoupling of impurities.- To dissect the RKKY 
regime, Jj > Jf, we take a central block of M spins which 
contains both impurities (see Fig. |TJb)) and compute its von 
Neumann entropy S(p(M)). Deep in the RKKY regime, 
Jj 3> Jf, the two impurities form a singlet and their quantum 
state becomes pure and decouple from the rest. This can be 
modeled by having two decoupled impurities in a singlet state 




FIG. 3: (Color online) (a) S{p %L l ) as a function of Ji. (b) Entan- 
glement between each impurity and its bulk vs. Ji. (c) Entanglement 
between each impurity and its opposite bulk vs. Ji. (d) Entangle- 
ment between the two bulks vs. J/. In all figures Nl, = Nr = 200. 
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FIG. 4: (Color online) The first derivatives with respect to Ji as a 
function of Ji for J' — 0.4 and two system sizes N = 240 and 
N = 400: The derivative of (a) Af(l L ,B L ); (b) Af(l L ,B R ); (c) 
M{Bl, Br), (d) The error parameter e vs. Ji for a system of N = 
240 and J' = 0.4. 



together with an effective system of length N — 2 formed by 
the left and right bulks as shown in Fig. [TJ C )- The effective 
interaction Hj* between the two bulks can be determined by 
a Schrieffer- Wolff transformation 12011 with the result that 



H 



eff 



J' 2 J X 

2Jj 



772 T2 
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2J/ 



2J/Ji 

(a 2 L .af + a*.a R - a£.a R - a 2 L . ( r 3 fl ),(4) 



where we have included also the modified boundary interac- 
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FIG. 5: (Color online) (a) Kondo length versus Ji. (b) Scaling of £k 
vs. 1/J' for different J/s. In both figures N L = N R = 200. 



Ji 


-3.00 


-2.50 


-2.00 


-1.50 


-1.00 


-0.50 


0.00 


a(J/) 


4.1939 


3.9976 


3.8861 


3.1598 


2.6330 


2.1262 


1.7753 



TABLE I: The values of a for different Ji's in a chain of Nl 
N r = 200. 



tion within each chain. Having effectively removed the im- 
purities, the central block can now be considered as that of 
M — 2 sites with density matrix p e g(M — 2). Given this ef- 
fective model one can compare the von Neumann entropy of 
the central block, i.e. S(p e s(M — 2)), with the one from the 
original system, i.e. S(p(M)), and only for Jj ^> Jf do we 
expect that the two entropies coincide for all M > 2. As a 
test we introduce an error parameter £ which quantifies the 
average difference between the von Neumann entropy of the 
central block for the original and effective chain, 

1 N 

W, Ji) = J^— 4 E \ S (P( M )) - S(Peff(Af - 2))|. (5) 

M=4 

In Fig.Htd) we plot £ as a function of Ji for J' = 0.4. Surpris- 
ingly, as revealed by the exponential decay of £, the impurities 
decouple from the bulk already for intermediate values of Jj. 

Kondo screening cloud.- We shall finally address the impor- 
tant issue of the size of the Kondo screening cloud. Follow- 
ing Ref. |@] we first trace out two identical blocks of spins, 
one from each bulk as shown in Fig. [TJd), and then compute 
the negativity between the impurities and the rest of the sys- 
tem. The negativity is found to decay exponentially with the 
number of spins which are traced out. We take the length be- 
yond which the negativity is less than a threshold value, here 
choosen as 0.01, to define the Kondo length 12111 . and plot it 
as a function of Ji in Fig. [2 a). As seen in the figure, £k in- 
creases as one tunes Jr from (with two independent single- 
impurity Kondo clouds) to large negative values (where the 
cloud is that of an exactly screened two-channel spin-1 Kondo 
model). Interestingly, increases with small positive values 
of J j and takes a small maximum at Jf. This indicates that 
the QPT is signaled also in the Kondo screening length at fi- 
nite system sizes. 

We have also investigated the scaling of the Kondo length 
£k as a function of the Kondo coupling J' for different val- 



ues of Ji < 0. In Fig. |5jb) £,k is plotted as a function of 
1/ J' in a semi-logarithmic scale. As £k approaches the total 
length of the system (N = 400) the curves for negative Jj 
start bending, indicating that the deviation from linearity is a 
finite-size effect. One thus infers from Fig. 0b) an exponen- 
tial dependence ~ e a l J in the large-volume limit, with 
a depending on Jj. Some values of a are shown in TABLE 
I. Extrapolation to large negative values of Jj suggests that 
a{— oo) « 2.5a(0) which is compatible with the results of 
the analysis carried in Ref. I23I1 . With the Kondo temperature 
T K related to £ K by T K = hv/k B £ K HI, our entanglement 
probe of the Kondo cloud thus shows an exponential suppres- 
sion of the Kondo temperature with the size of the effective 
impurity spin (i.e. spin-1 rather than spin- 1/2). This phe- 
nomenon, known as Kondo resonance narrowing l23tl . was 
implicitly touched upon in the original work on the TIKM (JI] . 
While the problem has been revisited recently in the context 
of the two-orbital Anderson model 12211 and for magnetic ions 
with large Hund's coupling I23I1 . ours is the first display of the 
effect extracted from quantum entanglement. 

Conclusion.- We have introduced a spin-chain model rep- 
resenting the two-impurity Kondo model and investigated 
its properties by applying two complementary entanglement 
measures borrowed from quantum information theory, nega- 
tivity and von Neumann entropy. This novel approach is con- 
ceptually simple and can easily be implemented numerically 
via a DMRG code. As we have shown in this Letter, it enables 
one to faithfully recover highly nontrivial features of the two- 
impurity Kondo model, including the existence of a quantum 
phase transition at a critical RKKY coupling [|2j]. Importantly, 
it makes possible, for the first time, a precise probe of how 
the elusive Kondo cloud depends on the Kondo and RKKY 
couplings, strikingly showing the effect of Kondo resonance 
narrowing within a controlled nonperturbative formalism. We 
expect that our approach can be exploited for generic quantum 
impurity problems and that it will prove increasingly useful 
with future applications. 
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